Breakup of small aggregates driven by turbulent hydro dynamic stress 
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Breakup of small solid aggregates in homogeneous and isotropic turbulence is studied theoretically 
and by using Direct Numerical Simulations at high Reynolds number, Re\ ~ 400. We show that 
turbulent fluctuations of the hydrodynamic stress along the aggregate trajectory play a key role in 
determining the aggregate mass distribution function. Differences between turbulent and laminar 
flows are discussed. A novel definition of the fragmentation rate is proposed in terms of the typical 
frequency at which the hydrodynamic stress becomes sufficiently high to cause breakup along each 
Lagrangian path. We also define an Eulerian proxy of the real fragmentation rate, based on the joint 
statistics of the stress and its time derivative, which should be easier to measure in any experimental 
set-up. Both our Eulerian and Lagrangian formulations define a clear procedure for the computation 
of the mass distribution function due to fragmentation. Contrary, previous estimates based only on 
single point statistics of the hydrodynamic stress exhibit some deficiencies. These are discussed by 
investigating the evolution of an ensemble of aggregates undergoing breakup and aggregation. 

PACS numbers: 47.27-i, 47.27.eb, 47.55.df 



Turbulence has a distinct influence on the aggrega- 
tion of colloidal and aerosol particles. It not only en- 
hances the rate of collision among particles, i.e. by in- 
ducing high velocity differences and preferential concen- 
tration within the particle field [TJ [2], but also it cre- 
ates hydrodynamic stress that can cause restructuring 
and breakup of aggregates [3 , a phenomena macroscop- 
ically expressed in shear thinning in dense suspensions 
[4]. Breakup of small aggregates due to hydrodynamic 
stress in turbulence is of high relevance to various appli- 
cations, e.g. processing of industrial colloids, nanoma- 
terials, wastewaters, and sedimentation of marine snow 
|31 El- In a mean field situation, aggregation-breakup 
dynamics are described by the Smoluchowski equation. 
Defining n^(t) = N^(t)/No where N%(t) is the number 
concentration of aggregates consisting of £ primary parti- 
cles, and Nq = / °° d^ £A^(£), the Smoluchowski equation 
reads as 
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where (j) = ^na^No is the solid volume fraction, a p is 
the radius of the primary particle (assumed monodis- 
perse and spherical), and fc^/ is the aggregation rate. 
Breakup is accounted for by the fragmentation rate f% 
Determining these functions is not easy and despite 
considerable efforts [6 a basic understanding of breakup 
dynamics is still lacking. A reason for this is the complex 
role of turbulence and the way it generates fluctuating 
stress to which an aggregate is exposed to. 



The main issue we investigate in this paper is how to 
define and measure the fragmentation rate fa in a turbu- 
lent flow. The outcomes of our analysis are manifold: i) a 
Lagrangian and an equivalent Eulerian definition of the 
fragmentation rate can be derived, that fall off to zero 
in the limit of small aggregate mass, while they have a 
power-law behavior for large masses; ii) the power-law 
tail description is crucial to obtain a steady-state aggre- 
gate mass distribution when considering the full aggre- 
gation/breakup dynamics; iii) turbulent fluctuations al- 
low for a broad asymptotic mass distribution, while a 
much narrower distribution of aggregates is obtained in 
the laminar case. 

We adopt the simplest possible framework [6|, and con- 
sider a dilute suspension of aggregates in a stationary 
homogeneous and isotropic turbulent flow. We con- 
sider very small aggregates - much smaller then the Kol- 
mogorov scale of the flow, in the range of 25 to 100 fim 
for typical turbulent flows, with negligible inertia. Fur- 
ther, we assume that aggregates concentration is such 
that they do not modify the flow. Hence, their evolution 
is identical to that of passive point-like particles (unless 
extreme deviations from a spherical shape are present, 
e.g. elongated fibers). Moreover, the aggregates are brit- 
tle and breakup occurs instantaneously once being sub- 
ject to a hydrodynamic stress that exceeds a critical value 
[7]- For small and inertialess aggregates, the hydrody- 
namic stress exerted by the flow is ^(e/^) 1 / 2 , where \i 
and v are the dynamic and kinematic viscosity, respec- 
tively and e is the local energy dissipation per unit mass. 
Thus, the key role is played by the turbulent velocity gra- 
dients across the aggregate which are known to possess 
strongly non-Gaussian, intermittent statistics [8]. 
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FIG. 1: Pictorial evolution of the energy dissipation e along 
an aggregate trajectory. Starting to record the stress at time 
to , the exit-time r and the diving-times T\ , T2 , . . . for a given 
threshold £ cr are shown. 



Let £ C r(0 be the critical energy dissipation needed to 
break an aggregate of mass £. In the simplest case of 
a laminar flow, e cr relates to the critical shear rate for 
breakup as G cr ~ {s CT /v) 1 ^ 2 . Earlier works [7]|9] support 
the existence of a constituent power-law relation for e cr 
implying that larger aggregates break at a lower stress 
than smaller ones: £ C r(0 = where the ex- 

ponent q is related to the aggregate structure and £ s is 
the characteristic aggregate mass. An equivalent rela- 
tion exists for small droplets breaking at critical capillary 
number, e CY ~ a 2 /(/ip£ 2//3 ) (here £ is the droplet volume, 
and a the interfacial energy). Hence, fa can be equally 
formulated in terms of a critical dissipation / £cr (£). 
In this contribution we propose to define the fragmenta- 
tion rate or f £cr in terms of a first exit-time statistics. 
This amounts to measure the fragmentation rate using 
the distribution of the time necessary to observe the first 
occurrence of a local hydrodynamic stress strong enough 
to break the aggregate. An operational formula of f £cr 
reads as follows: (i) seed homogeneously a turbulent, sta- 
tionary flow with a given number of aggregates of mass 
£; (ii) neglect those aggregates in regions where the hy- 
drodynamic stress is too high (e > e cr ); (hi) from an 
initial time to, selected at random, follow the trajectory 
of each remaining aggregate until it breaks, and count 
the total number of breaking events in a given time in- 
terval [r, r + dr]. The time r is the first exit-time: for 
an aggregate initally in a region with e < s cr , r is the 
time it takes to the hydrodynamic stress seen by the ag- 
gregate along its motion to cross the critical value e cr at 
first opportunity, Fig. [I] The fragmentation rate is the 
inverse of the mean exit- time: 
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where V £cv (r) is the distribution of first exit-time for a 
threshold s cr , and (-) ex is an average over V £cT (r). The 
definition in Eq. Q is certainly correct but difficult to 
implement experimentally as it is needed to follow aggre- 
gate trajectories and record the local energy dissipation, 
something still at the frontier of nowadays experimental 
facilities [10]. The question thus arising is if we can ob- 



tain a proxy for the fragmentation rate which is easier to 
measure. Given a threshold for the hydrodynamic stress 
£ cr , one can measure the series T±, T2, ... of diving-times, 
namely the time lags between two consecutive events of 
istantaneous stress crossing the threshold s cr along the 
aggregate motion, Fig. [I] In [IT], it was proposed to es- 
timate the fragmentation rate as the inverse of the mean 
diving-time, = l/(T(e cr )). An important result is 
that (T(e cr )) can be obtained using the Rice theorem for 
the mean number of crossing events per unit time of a 
differentiable stochastic process across a threshold [12]. 
Hence, 
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1 _ / °°rig ep 2 (e cr ,e) 



(3) 



Here the numerator is the Rice formula giving the mean 
number of crossings of e cr in terms of the joint proba- 
bility of dissipation and its time derivative p 2 (£,£); the 
denominator is the measure of the total time spent in 
the region with e < e cv . Notice that the integration in 
the numerator goes only on positive values in order to 
consider only up-crossing of the threshold e cr [12 . An 
obvious advantage of Eq. ([3| is that it is quasi- Eulerian: 
it does not require to follow trajectories, since it depends 
only the spatial distribution of dissipation and of its first 
time derivative in the flow. Expressions ^ and (|3| are 
not strictly equivalent. A direct calculation of the mean 
exit-time in terms of the distribution of diving-times gives 
indeed (r(e cr )) ex = {T 2 (e cr )) /[2(T(e cr )>], which relates 
the mean exit-time to moments of the diving-time. 
Once a definition from first principles is set-up, we pro- 
ceed to measure the fragmentation rate for aggregates 
convected as passive point particles in a statistically ho- 
mogeneous and isotropic turbulent flow, at Reynolds 
number Re\ ~ 400. Details on the Direct Numeri- 
cal Simulations (DNS) of Navier-Stokes equations with 
2048 3 grid points and Lagrangian particles are in [2 . The 
present analysis is obtained averaging over 6 x 10 5 trajec- 
tories, recorded every 0.05^, where is the Kolmogorov 
time of the flow. 

Figure [2] shows the fragmentation rate measured from 
the DNS data following the evolution of the velocity gra- 
dients along aggregate trajectories. The exit-time mea- 
surement Q and its Eulerian proxy ([3| show a remark- 
able, non trivial behavior for small values of the critical 
threshold, i.e. for large aggregate mass. In this region, 
there is a competing effect between the easiness in break- 
ing a large aggregate and the difficulty to observe a large 
aggregate existing in a region of low energy dissipation. 
As a result, the estimated fragmentation rate develops a 
quasi power-law behavior for small thresholds. On the 
other hand, for large thresholds the super exponential 
fall off is expected. It is the realm of very small aggre- 
gates that are broken only by large energy dissipation 
bursts. The exit- time Q and diving time Q measure- 
ments are very close and we therefore consider the latter 



3 




FIG. 2: The normalized fragmentation rate fs cr Tri versus the 
normalized energy dissipation e C r/{s). Filled circles are the 
definition |2]), measured up to thresholds where statistical 
convergence of exit-times is obtained; the solid line is f^ cr . 
Squares are // cr , while crosses are // c J r . (Bottom inset) Joint 
distribution p2(e,i). The continuous line is the dimensional 
estimate e ~ e/r v (e). (Top inset) Numerator of Eqs. (|3| - 
(pj; curves colors are the same of the main figure. 




FIG. 3: Time evolution of Io for df = 2.4 and q = 0.36 with 
breakup rates ff cr (squares); // cr (circles); // c J r (triangles). 
Runs with different solid volume fractions (j) are shifted. Inset: 
Io at steady state as a function of (j> (same symbols). At large 
0, the model with fragmentation (13]) relaxes to the predicted 
scaling curve Io ~ ( j ) 1 /0-+x/<i-z/djT^ given by the solid line. 

Here, £ s = 10 4 implying that a/a p ~ ^J df ~ 50 is the value 
of the characteristic aggregate size. 



a very good proxy of the former, the real fragmentation 
rate. The main advantage of the estimate Q is the very 
high statistical confidence that can be obtained since it 
is a quasi-Eulerian quantity. Moreover it gives a reliable 
estimate also in the region of large thresholds where the 
convergence of exit-time statistics, requiring very long 
aggregate trajectories, is difficult to obtain. 
Starting from Eq. simple models can be proposed for 
the statistics of e, whose direct measure requires a very 
high sampling frequency along Lagrangian paths. First, 
as it appears from the bottom inset of Fig. |2j the dissi- 
pation and its time derivate are significantly correlated. 
Scaling on dimensional grounds suggests e ~ e/r^e) 
where r^e) ~ (v/e) 1 ^ 2 is the local Kolmogorov time. 
It follows that the joint PDF p2{s 1 e) can be estimated 
a 1 sp2(s 1 e) = \p{e)8{\e\— e/r^e)) where p(s) is the prob- 
ability density of energy dissipation. Prefactor 1/2 ap- 
pears since for a stationary process e is positive or neg- 
ative with equal probability. Plugging this expression in 
Eq. ^ gives 

j, j _ 2^crP(^cr) / T r](^cr) 

h «~ J*"dep(e) ' U 

We refer to it as Closure I. A different approach was 
proposed in [6] . It assumes that active regions in the flow 
where e > e cr engulf the aggregates at a rate ~ l/r v (e), 
which results in 

f™dep(e)/T v (e) 
h " Jtdep(e) ■ U 

We refer to it as Closure II. Both models share the ad- 
vantage of being fully Eulerian and based on the spatial 



distribution of the energy dissipation only. In Figure [2j 
the fragmentation rates obtained from the closures I and 
II are also shown. Both of them reproduce the correct 
behavior for large values of the critical dissipation but 
deviate for small ones. The reason for such discrepancy 
is made clear in the top inset of Fig. |2j where the numer- 
ator of Eqs. Q and ([5| are shown. It appears that, for 
small values of the critical stress £ cr , // underestimates 
the number of breakup events, while the numerator of 
fl 1 saturates to a constant value. 

J t C r 

Since in experiments breakup a fortiori takes place to- 
gether with aggregate recombination, we explore how the 
actual fragmentation rate Eq. Q and the two closures, 
Eqs. Q-^, influence the time evolution of an ensem- 
ble of aggregates n^(t). At this purpose, the Smolu- 
chowski equation (fil), subject to the initial condition 
n^(0) = 5(£ — 1), is evolved in time. To model ag- 
gregation we use the classical Saffman- Turner expres- 
sion % e = A)Av?(£ 1/d/ + f 1/df ) 3 , where D is an 
O(l) constant and df is the fractal dimension that re- 
lates the collision radius of an aggregate to its mass, 
a/a p = Breakup is assumed to be binary and 

symmetric, g^^ = 25({; — £'/2), which despite its sim- 
pleness represents well the quality of the evolution. To 
quantify our findings we consider the second moment ob- 
servable readily accessible from static light scattering [7] . 
Fig. [3] shows the time evolution of Io for typical values 
of df and q found in turbulent aggregation of colloids 
[3]. After an initial growth period, curves obtained with 
Eqs. Q and ^ both relax to a steady state. At small 
solid volume fraction (/>, the two models nearly overlap, 
whereas at larger <\> Closure II underestimates Io. This 
behavior is confirmed in the inset of Fig. [3] that shows Jq 
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at steady state for both models and for different values of 
(j). Clearly, at large <fi other phenomena, i. e. modulation 
of the flow due to the particles, may occur which, how- 
ever, is beyond the scope of the present work. Closure 
I shows a very different behavior. At small an evolu- 
tion similar to previous cases is observed. However, by 
increasing <j> a drastic change appears and Iq diverges, a 
direct consequence of the presence of a maximum in the 
shape of // cr , see Fig. 2j 

The time evolution ofthe number concentration 
governed by Eq. ([3|, is further examined in Fig. [4j Here, 
we compare the steady distribution obtained in the tur- 
bulent flow with that of a laminar flow for different values 
of the solid volume fraction. In the turbulent case, n^(t) 
rapidly grows and reaches a stationary state whose peak 
mode is controlled by the magnitude of aggregation, i.e. 
the solid volume fraction (f). Increasing <ft causes the mode 
to broaden and to shift to the right as aggregation gets 
more pronounced. The distribution thus gradually moves 
into the region where f £cr assumes power-law behavior. 
A numerical fit of the left fragmentation tail in Fig. [2] 
gives ft ~ fx/*, with x = 0.42 ± 0.02 (dashed curve m 
Fig. [4| . Using this latter expression in Eq. ([!]) , one can 
derive a scaling relation for integral quantities of n^(t) at 
steady state [14]. Such scaling is reported in the inset of 
Fig. [3] On the other hand, in the laminar case where a 
uniform shear rate governs the breakup, the steady state 
distribution is much narrower and shows multiple reso- 
nant modes. These are due to the sharp onset of breakup 
once the aggregates grow larger than the characteristic 
aggregate mass £ s . These results clearly demonstrate the 
strong influence of turbulent fluctuations on the statisti- 
cally stationary mass distribution function. 
We have presented a study of the fragmentation rate of 
small and diluted aggregates in turbulent flows at high 
Reynolds number. We have introduced a novel expression 
for the fragmentation rate in terms of the exit-time statis- 
tics, which is a natural way of measuring first-order rate 
process. Also, a purely Eulerian proxy based on Eq. (|3| 
provides a very good approximation to the actual frag- 
mentation rate measured from our DNS. Remarkably, a 
steady state in the full breakup-aggregation process is 
crucially determined by the left tail of the fragmentation 
rate, i.e. by events of low energy dissipation. Our in- 
vestigation puts the basis for many developments, such 
as the stability of the Smoluchowski evolution using the 
measured fragmentation rates in experiments, and the 
extension to the case of inertial aggregates. In such case, 
the correlation between the hydrodynamic shear and the 
Stokes drag may result in a non-trivial breakup rate de- 
pendency on the degree of inertia. Future work aims 
to introduce spatial fluctuations in the mass distribution 
caused by local breakup, a research path still poorly ex- 
plored. 
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FIG. 4: Stationary mass distribution n^(t) with df = 2.4 and 
q = 0.36, for different solid volume fractions (j). Distributions 
obtained in the turbulent flow at Re\ ~ 400 (a) are compared 
to those of the laminar cases (b) , with the uniform shear rate 
£ l cr m — (e). In (a), we plot (dashed line) as given by Eq. ([3]), 
assuming a power-law behavior in the limit of large aggregate 
mass. 
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